NASA Technical Memorandum 106805 



Microgravity Isolation System Design: A Modem 
Control Synthesis Framework 


R.D. Hampton 
McNeese State University 
Lake Charles, Louisiana .. . 

C.R. Knospe and RE. Allaire 
University of Virginia 
Charlottesville, Virginia 


C.M. Grodsinsky 
Lewis Research Center 
Cle veland, Ohio 


■ mm 1 w i a - r uaE it. 1 ^.g p awMBBaaBBgT!*raBt!h ■ 1 ~^~ ; 


December 1994 



National Aeronautics and 
Space Administration 


(NAS A- TM- 106805) MICROGRAVITY 
ISOLATION SYSTEM DESIGN: A M00ERN 
CONTROL SYNTHESIS FRAMEWORK (NASA. 
Lewis Research Center) 34 p 


N95-18197 


Unci as 


G3/01 0034990 



Microgravity Isolation System Design: 

A Modern Control Synthesis Framework 

R. D. Hampton , 1 C.R. Knospe , 2 P.E. Allaire , 2 and C.M. Grodsinsky 3 


Abstract 

Manned orbiters will require active vibration isolation for acceleration-sensitive 
microgravity science experiments. Since umbilicals are highly desirable or even indispensable for 
many experiments, and since their presence greatly affects the complexity of the isolation 
problem, they should be considered in control synthesis. In this paper a general framework is 
presented for applying extended H 2 synthesis methods to the three-dimensional microgravity 
isolation problem. The methodology integrates control- and state frequency weighting and input- 
and output disturbance accommodation techniques into the basic H 2 synthesis approach. The 
various system models needed for design and analysis are also presented. The paper concludes 
with a discussion of a general design philosophy for the microgravity vibration isolation problem. 

Introduction 

Although many scientists have planned or conducted materials processes and fluid physics 
science experiments designed for a weighdess environment, the cunendy available facilities have 
proved far from ideal. Evacuated drop towers can provide only a few seconds of 
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"weightlessness" at levels on the order of 10* go (where go is the gravitational acceleration at sea 
level). Aircraft flying low-gravity parabolic trajectories can extend the time to about 15-20 
seconds, and sounding rockets can provide several minutes of a microgravity environment; but the 
goal of providing days, or even hours, for microgravity research has proved elusive. It was once 
hoped that the Space Shuttle could provide the desired environment; but such factors as manned 
activity, machine and structural vibrations, and thruster firings for orientation or reboost have 
resulted in acceleration levels generally unsatisfactory for the designed experiments. (Background 
excitations have been measured in the milli-g range.) In fact, the data from many experiments 
have been found unacceptable due to the poor acceleration environment 

Due to the low frequencies of greatest concern (below about 10 Hz) the isolation problem 
is a largely unfamiliar one to vibration engineers; the requirement of a comer frequency of about 
10-3 Hz is particularly vexing. Passive isolation is incapable of solving the isolation problem in 
this region; and even should a sufficiently soft spring be physically realizable, it could not isolate 
against direct disturbances to the payload. If the payload is tethered (e.g., for evacuation, power 
transmission, cooling, or material transport), a passive isolator cannot provide isolation below the 

comer frequency imposed by the umbilical stiffness. 

An active isolator (such as a magnetic suspension system) that merely possesses a low 
positive stiffness fares no better in the presence of an umbilical, for the same reasons. And if the 
control system seeks to lower the comer frequency by adding negative stiffness, to counteract the 
umbilical's stiffness, the system will (at best) possess almost no stability robustness. In the face of 
the usual umbilical nonlinearities and uncertainties, this situation is clearly unacceptable. At very 
low frequencies the rattlespace constraints become limiting [1,2], so that any isolation system 
must have unit transmissibility in that region. In short, a microgravity isolator must be active, and 
it must be capable of dealing with the particular frequency-dependent complexities accompanying 

a tethered payload and a restrictive rattlespace. 

The available acceleration data clearly point to a need for three-dimensional isolation [3]. 
Classical control design methods are not well-suited for handling such problems; modem control 
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methods provide a much more natural setting, opening up to the designer the power of the 
developing robust control synthesis and -analysis tools, along with a variety of well-tested and 

progressive computational software packages. 

The well-known H 2 synthesis [i.e., LQR ("Unear Quadratic Regulator") or LQG("Linear 
Quadratic Gaussian")] methodology is one such modem control method. It can readily provide an 
optimal feedback controller for a linearized plant (i.e., payload plus umbilical) subject either to no 
exogenous input (LQR case) or to white noise disturbances only (LQG case). An optimal control 
found by H 2 synthesis minimizes a quadratic ("energy-type") cost function, or performance index. 
Such a performance index is quite appropriate for the microgravity isolation problem, since it 
allows penalizing both the control energy required for isolation and the vibrational energy of the 
payload. Unfortunately, however, the application of this synthesis method to practical problems 
has been plagued by robustness difficulties. Granted, the standard LQR solution provides 
excellent robustness guarantees for the single-input-single-output (SISO) problem ([4], pp. 70-74) 
and also yields guarantees (though less useful) for the multiple-input-multiple-output (MIMO) 
problem [5]. But the addition of a state observer to the controller (as is usually neccessary for 
practical problems) removes these robustness guarantees [6]. This fundamental practical concern 

has led to a common skepticism regarding H 2 synthesis. 

There exist extensions to H 2 synthesis, however, which can resolve the robustness issues. 
The disturbance-accommodation and frequency-weighting techniques contributed, respectively, 
by C. D. Johnson ([7], 1968; [8], [9], 1970; [10], 1971) and N. K. Gupta ([1 1], 1980) have 
proved to be highly useful in this regard. In fact, they provide the fundamental additional tools 
needed for solving practical controller design problems. These two extensions lead to augmented 
state equations which still allow for problem solution by the familiar Hj synthesis machinery. 
Recent investigations have examined the effect of the frequency-weighting extension on system 
robustness [12,13], and the dual relationship between frequency weighting and disturbance- 
accommodation [14], Additional extensions have also been proposed [15]. 
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The utility of "extended H 2 synthesis” for the tethered microgravity vibration isolation 
problem has been clearly demonstrated by recent studies [16.171. Extended H 2 synthesis has 
been used effectively to develop SISO and SIMO (single-input-multiple-output) controllers for a 
realistic one-dimensional microgravity vibration isolation problem, using a smart form of 
acceleration feedback. The resulting closed-loop system exhibited excellent stability- and 
performance robustness guarantees, including a high degree of robustness to umbilical parametric 

uncertainty. 

The present paper will give a general framework for controller design by the extended H 2 
synthesis method, for the microgravity vibration isolation problem. Following a summary of the 
basic H 2 synthesis approach, the paper will describe how to incorporate the control- and state 
frequency weighting and input- and output disturbance accommodation extensions into the 
synthesis procedure. Control noise will also be included. General guidelines will be presented for 
effectively integrating these extensions into the design procedure. 

The development below is specifically tailored to the microgravity isolation problem, but 
the mathematics are fully applicable to any problem that has the appropriate (very general) 
mathematical description. Four fundamental system models will be presented to aid the designer 
in visualizing the design effort. Only the synthesis procedure will be detailed here; the analysis 
techniques used for controller evaluation will be detailed in a later work. 


Basic H 2 Synthesis Review 

A generic microgravity vibration isolation system is depicted below in Fig. 1. A payload, 
such as a microgravity science experiment, is acted upon by actuators (typically non-contacting) 
that are commanded by a control system. This control system uses measurements, such as 
payload positions and accelerations, to develop the control signals, typically currents or voltages. 
The objective of H 2 synthesis is to find a control signal that minimizes the weighted sum of the 
two-norm of the control energy and of the states, subject to the linearized system equations of 
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motion. This control signal will be found to be dependent only on the past accumulative 
measurement information, for a system excited only by zero-mean white Gaussian noise. 


Actuator 


Payload 


Sensor 


A 




Control signals 



Measurements 


Figure 1.— Vibration isolation system. 


Specifically (and a bit more mathematically), to use H 2 synthesis the system equations of 
motion must first be linearized and expressed in the following (standard) state space form. 


x = Ax + Bu + E s w t 
y = Cx + Du 
l = y+E n w n 


(la) 

(lb) 

(lc) 


i is the state vector, y is the output vector, z is the measurement vector, u is the control vector, E s 
and E n are selection matrices, and and jy, are process- and sensor-noise vectors, respectively. 
For the microgravity vibration isolation problem, the process noise (j^) models Ac disturbances 
acting on the payload, either directly (c.g., air currents, fluid flow, or experiment-mounted 
rotating machinery) or indirectly (i.e., through the umbilicals). The sensor noise GO models Ac 
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electrical or mechanical noise that contaminates the state measurements. In general, not all states 
will be measurable (i.e., rank C dim i}. 

Let the initial conditions on the state vector be *(0) = iQ (although these intitial conditions 
will not appear in the final control solution); let ais > an ^ &n ^ independent and bounded (as is 
reasonable, since unbounded states and infinitely large noise are not physically possible), and let 
iQ be Gaussian, and a' n and njy, zero-mean white Gaussian, for technical reasons. The power of 
the process- and sensor-noise vectors can be expressed mathematically by 

covLa^O^Ct )l = V r ; 5 (r - x ) and cov[^(r), xjx )} = Vj5 (t-% ) (2) 

([18], p. 272). Assume that [A#] and [A,E, V/' 2 ) are stabilizable, where V t = V, lli V," 2 * (the 
asterisk here means "conjugate transpose"); and that [CM is detectable ([19], p. 226). These 
requirements mean, respectively, that a stabilizing controller exists, and that the available 
measurements are sufficient for its implementation. Let V 2 and V 3 be positive semidefinitc (PSD) 
and positive definite (PD), respectively for reasons of solution existence. That is, there need not 
be any process noise, but there must be at least some noise in all measurement channels (as there 

will be) if an optimal control solution is to exist 

The H 2 synthesis design method uses a quadratic performance index. 


2T 


/ J 

1 

J* 

1 

X 

\1 

•e([i r /] 

k ' 



l W 2 W > J 

u 

I/I 


dt 


(3) 


where W t , W 2 , and W 3 are weighting matrices. These weighting matrices, assigned by the 
designer, allow him to place a relative importance on the two- norm of each state (using W t ) and 
of the control (using W 3 ). W 2 allows him to assign cross weightings. (These cross-weightings are 
not generally used for the basic H 2 synthesis problem, but they become important with some of 
the extensions.) W t is PSD and W 3 is PD ([18]; pp. 272, 276), again for reasons of solution 
existence. "5" is the expected-value operator, needed since the system is excited stochastically by 
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Uj,. The cost rate functional form for J (with " lim — " ) is used to allow for the white noise 

disturbance a; Otherwise the performance index would be infinite. 

Let an admissible control u(0 be one that depends only on the past accumulative 

observation data. That is, u(t) has the form 

{£(/) = u[f,Z(/)] < 4a > 

where Z(t) = {z(x),0£x£f}. 

(For more general conditions on admissibility, see [18], p. 272.) The objective is to find an 
admissible control function u (t) which minimizes the cost / with respect to the set of admissible 

control functions u(t) subject to the dynamic constraint (la,b,c). That is, the optimal control 
solution must exist and be realizable, must minimize the cost functional specified by the designer, 
and must take into consideration the system equations of motion. 

The solution is well-known, and is summarized as follows: 
m *(0 = -* x (0 

where I is an estimate of x using a Luenberger observer 
([18], pp. 288-289) having observer gain matrix L 


k=w; 1 {b t p+wJ), 

(5b) 

P is the unique PD solution 


to the Algebraic Riccati Equation (ARE) 


PA + A T P-(PB + W 1 )Wf'(PB + W 2 ) T +W, =0, 

(5c) 

l=qc t (e,v j e t J', 

(5d) 

Q is the unique PD solution to the ARE 

AQ + QA T -QC T {E n V 3 El)~ 1 CQ + E,V 1 E] =O f 

(5e) 


P exists if [AR] is stabilizable and [CA) is detectable 
or if the system is asymptotically stable, and 
Q exists if [A, E,Vj m ) is stabilizable and {CA } is detectable 
or if the system is asymptotically stable. 


7 



Note that this control signal is developed by simply applying constant (negative) feedback gains K 
to estimates of the system states. These estimates are themselves optimal in that they are the 
"closest" to the actual states, in terms of the expected value of the rms of the estimate error. They 
are produced in the controller from the control signal and the measurement vector, using constant 
observer gains L. Fig. 2 below presents this standard optimal controller (e.g., cf. [20], pp. 356, 



Figure 2. — Block diagram of system with ^-optimal controller. 


Extensions to H2 Synthesis 

Frequency Weighting 

For the active microgravity vibration isolation problem, payload accelerations are of much 
greater concern at some frequencies than at others. Accelerations at higher frequencies can be 
handled passively; and very low-frequency accelerations correspond to such large displacements 
that they are essentially unisolable, due to practical rattlespacc constraints. Rattlespace 
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constraints also require that the relative displacements between space platform and payload be 
kept to a minimum at low frequencies. Control is needed at lower frequencies, where the plant is 
best-known and where the major isolation effort is desired. At higher frequencies, however, 
excessive control can excite unmodeled higher modes of the plant Consequently it is desirable, 
through the performance index, to be able to penalize control strategies in a frequency-dependent 
fashion. This can be achieved by weighting the states i and the control u in the cost rate 
functional so that the weightings are frequency-dependent The latter, extended, H 2 synthesis 
problem will be seen to have the same form as the former one. It will have the simple difference 
dial the various system matrices now will be "augmented" to take into account the additional 

"pseudo" states required by the frequency-weighting extension. 

Let * be considered to be the input to a filter #s) of which * is the output, and let 


tfis) have a state-space representation defined by [A,, B lt C,, D { ). 
That is, fts) = CjisI-AjY'B'+D,]. 

Then 1/ = A t z 2 + B t x 

x = C 1 z i +BjX 


( 6 ) 

(7a) 

(7b) 


expresses X in terms of j. employing pseudostates z t . Similarly, if U is considered to be the input 
to a filter *s) of which u is the output, and if ^s) has a state-space representation defined by 


[A 2 , B 2 , C 2 , D 2 ),u can be expressed in terms of u, employing pseudostates 12 : 


z 2 — A 2 z L 2 B 2 U 
u = C 2 z 2 + D 2 u 


(8a) 

(8b) 


These frequency-weighted states (*) and controls (u) are now weighted (Le., penalized) by 


constant weighting matrices W t and W 3 , respectively. The resulting state equations and 

performance index are as follows: 
l x= l A l x+ *Bu+ 

y= l C I x + Du 
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(9a) 

(9b) 



(9c) 


where 


z = y+E n w 


‘J 



&F<“ r ([' *v{ 


% v 2 

V/ 



(9d) 



(9e) 


'A = 


A O O 
B, A, O 
0 O A 2 


(9f) 


l B = 



(9g) 


! C = [C O O ] 




DjW,D, 

CjW,D, 


O 


dJw i c, o 
c]w t c l o 
0 c]w,c 2 


(9h) 


(9j) 


(9k) 



O 

O 

cJw 3 d 2 


(9m) 


10 


II ! 


1 w 3 =[dJw 3 d 2 ] (9n) 

The optimal control u(t) wiU now minimize the weighted sum of the two- norm of \he frequency - 
weighted control energy and states. 


Input Disturbance Accommodation 

In the basic H 2 problem it was assumed that the process noise w, (i.e., the disturbance 
acting on the payload, whether directly, or indirectly via the umbilicals) is zero-mean white 
Gaussian. This, of course, will not generally be the case; the process noise will have some 
(known or unknown) non-white power spectrum. Let the process noise be modeled as/j, where 
jr is a stochastic disturbance with power spectral density S^(g>) = Sj (yo>) Sj ^jca). Defining 
H by $ //2 0'co) y 112 , one can consider £ to be the output of a filter excited by zero- 

mean white Gaussian noise u?, with power (i.e., cov^fO,^ )] = v i& ( r * x ))• 

In state-space form. 


= + (10a> 

/,=C,| J + D,w, (10b) 

such that H f(s) = C,[sl - A,)~' B, + D, <10c> 

Incorporating these new pseudostates (£ ,) into the state equations and petfonnance index yields 

the further augmented H 2 synthesis problem given below; 

2 x = 2 A 2 x + 2 Bu+ 2 E J w J 

(lib) 
(11c) 


y = 2 C 2 x + Du 


z = y+£„w„ 


2 J 


= £ 


lim — 
t-*~ 2T 




(lid) 
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where 


J'JEI 




kj 




'A 

0 

0 £,C,] 

B, 

A, 

0 

0 

0 

0 

‘ Aj 

0 

0 

0 

0 

A, _ 


\B 



O 


(He) 


(Ilf) 


dig) 


2 C = [C 000 ] 

r £ Al 



B, J 


! IV,= 


d]w,d, 

CjW,D, 

0 

o 


DjW,C, 

C]w,c, 

o 

o 


o o 
o o 
c]w 3 c 2 o 
o o 



o 

o 


c]w 3 d 2 



2 w 3 =[d t 2 w 3 d 2 ] 


(Uh) 


(llj) 


(Ilk) 


(11m) 


(tin) 


12 


The optima] -control solution to this problem will minimize the frequency-weighted cost functional 
as before, with the plant now considered to be subject to the specified colored noise disturbance. 

In actual space applications the power spectrum of the process noise may not be known. 
Fortunately, orbiter spectral vibration information need not be available for disturbance 
accommodation to be used. The disturbances can be assumed to have whatever form the designer 
finds useful. For example, if he desires the controller not to respond to process noise above some 
frequency range, he might choose to model Hj(s) as a lowpass filter. 


Output Disturbance Accommodation 

The same procedure can be employed to incorporate colored sensor noise into the 
extended H 2 synthesis problem. No sensor will have the white-noise contamination assumed by 
the basic H 2 problem. And the designer may even find it useful to shape the sensor noise filter in 
some non-physical way. For example, if payload acceleration measurements are known to be 
more accurate in one frequency range, and relative position measurements in another, he might 
choose his sensor noise filters appropriately to "inform" the observer of these facts. The resulting 
observer would tend to rely more heavily on the more accurate measurement(s) in a particular 


frequency range, in its state-reconstruction process. 

Let the sensor noise vector be£,« with selection matrix £,, where is a stochastically 
modeled disturbance with power spectral density S m {Jn) = sj!\jcd) . As with the mput 


disturbance/,, £, can be considered to be the output of a filter H„(s) excited by zero-mean white 
Gaussian noise Hi, with power V 3 (i.e., cov[ (f)» *!« (*)] = V^5(r — x)). In state-space form. 


i 2 =A n + 

/„ =C n $ 2 + D„w n 

such that H n (s) = C„ {si - A n )~ l B H + D„ 


(12a) 

(12b) 

(12c) 


For the extended H 2 synthesis problem with state- and control frequency weighting 


(pseudostates i ; and Zz* respectively), and with input- and output disturbance accommodation 
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(pscudostates £ , and £ 2 , respectively), the augmented state equations and the performance 

index are as follows: 

4 ±= 4 A 4 x + 4 Bu+ 4 E j 4 w j 

z=C 4 x+ 4 Du¥ 4 E„ 4 w h 


(13a) 

(13b) 



(13c) 


x 

h 


where 



I —2 




5 , 

5, 


4 a = 


A O 

B, A t 
O O 
O 0 
O 0 


O E,C, O 

OOO 

a 2 o o 

O A, O 

O O A„ 


4 B = 


B 

O 

B 2 

O 

O 


4 C = [C OOO E.C.] 
4 d = d 


(13d) 


(13c) 


(13f) 


(13g) 

(13h) 
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4 e = 


E,D, O 
O O 
0 O 
B. O 


O B, 


* J 



(13j) 


(13k) 




(13m) 

(13n) 



DJW 1 D 1 
cj Wj D, 

O 

o 

o 


Djw t C, o 0 0 

cj w l c 1 o 0 0 

o cjw 3 c 2 o o 

o 0 0 0 

O 0 0 0 


(13p) 


o 

0 



o 

o 


(13q) 



(13r) 



O' 

Vs 


is the autocorrelation matrix for 4 w s , 


(13s) 



the cross-correlationmatrix between 



andvv,, (13t) 
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and 4 Vj =V 3 ls\hc autocorrelation matrix for 


(13u) 


Control Noise 

The active microgravity vibration isolation system must perform well even when the actual 
system dynamics are less than perfectly described by the system model. One way to improve the 
isolation system’s robustness to parameter changes at the control inputs (i.e., to uncertainties in 
the actuator or plant model) is to add a process noise input (^) to the model’s control signal 
([21], p. 1-48). Recall that the H 2 synthesis machinery seeks to minimize the rms of the 
observation error. It does this by finding an observer gain matrix L that will optimally trade off 
the measurement uncertainties against the plant model uncertainties, in the state reconstruction 
process. Control noise will reduce the observer’s "confidence" in the plant model, so that the 
observer will "trust” its measurement data more and its plant model less. The resulting gain 
matrix L will sacrifice a degree of observation quality for improved observer robustness to plant 


model inaccuracies. The controller gain matrix K will be unaffected. 

Under these circumstances the state equations of motion, unaugmented by frequency 

weighting or disturbance accommodation, become 

x = Ax + £(« + vyJ + £,/, 

y = Cx + £>(u + vyJ (14b) 

r. , (14C) 

i =y+ E J K 


where £ and £ can be represented by filters in state-space form with white-noise inputs, as noted 

previously. Assume no cross-correlation between and or between and and let 

cov [^(r), H^(t)] = V r c 5(r-x). 


Using now the pre-superscript designator ”5” to indicate the appropriate state-space 
augmentation, the system equations change as follows: 

5 x = 4 x 

S A= 4 A, S B=B, 5 C=C, 5 D= 4 D, 


(16a) 

(16b) 
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*w t = 4 w t , s w 2 =w 2 , 5 W 3 = 4 Wj 


(16c) 



B 

O 

b 2 

O 
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E,D, 

O 

o 

B, 

O 


O 

O 

O 

o 


(16d) 


5 Wj =<i 


w c 

w. 


(16c) 



(160 


5 E„=[D E n D n ] 


( 16 g) 


5 


*5» 


y e o o 

O Vj o 

o o v 3 


(16h) 


5 


^2 = 


K 

o 

o 


o 

y 2 

Vj 


(16i) 


5 


V,= 


X 

o 


O' 

V, 


( 16 j) 


The basic tools are now in place for practical micro gravity vibration isolation system design, by 


extended H 2 synthesis. 


System Modeling 
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The H 2 synthesis problem is actually a two-fold design problem; the designer must 
determine a regulator gain matrix K, and also an observer gain matrix L, which together are used 
to comprise the optimal controller. The full augmented state vector must be used for the 
regulator sub-design problem; but pseudostates i, and i 2 (which occur due to frequency 
weighting) need not be reconstructed by the observer. They can simply be developed by passing 
the reconstructed state vector i and the control vector if, respectively, through the appropriate 
frequency-weighting filters. Consequently the observer sub-design problem can (though need 
not) be one of reduced older. It is helpful, then, to have different mathematical models for 
c once p tualizin g these two sub-ptoblems. These models will also differ, in general, from the basic 
plant model, which depicts the "actual" linearised plant <i.e.. the system without the controller) in 
state space form. This model will not include the frequency-weighting and disturbance- 
accommodation augmentations of the former. In addition to these three conceptual models of the 
system, them is a fourth model, which more properly falls under the category of analysis but 
should be kept in mind during the synthesis procedure. This "nominal analysis model" depicts the 
linearized and unaugmented plant with the synthesized controller attached. It is used, with 
various modifications, to analyze closed-loop system performance. 


Basic Plant Model 


The basic plant model (shown schematically in Fig. 3) simply presents the linearized 
differential equations of motion in a state space form useful to the H 2 synthesis machinery. Such a 
representation is given below. 


x = Ax + Bu + E,f 
y = Cx + Du 

z = y+E*f m 


(17a) 

(17b) 

(17c) 


No performance index is needed at this stage, since it is the frequency-weighted states and control 
which will be weighted relative to each other for the actual controller synthesis. For the 
microgravity vibration isolation problem a useful choice for the state vector would include relative 
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displacements, relative velocities, and accelerations. Weighting an acceleration more heavily in a 
frequency range would correspond roughly to a demand to increase the associated effective mass 
(or inertia) of the system . A similar correspondence can be drawn between relative-displacement 
weighting and the effective relative stiffness, and between relative-velocity weighting and the 
effective relative damping. 



Regulator Synthesis Model 

The regulator synthesis model adds frequency weighting and disturbance accommodation 
weighting filters to the basic plant model, and is the model actually used by the extended H 2 
synthesis machinery in designing the regulator. Fig. 4 portrays this model in block diagram form. 
Note that it has a similar form to the basic plant model, but that now augmented A, B, C, E S E H , 
W 2 and W 3 matrices are used, as indicated by the pre-superscripts. These matrices were 
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defined previously (16b through 16g). Also note that white noise vector u, replaces^, since the 
process noise power spectral information is now contained in matrices 5 A and 5 E t . 

The regulator synthesis model shows the system as viewed by the extended H 2 machinery 
in determining the matrix 5 P. This matrix is the unique positive definite solution to the following 

ARE, 

5 p 5 a+ 5 a t 5 p-( s p s b+ 5 w 2 ) 5 w 3 - 1 ( 5 p 5 b+ 5 w 2 ) t + 5 w,=o (18) 

J/> is used to find the regulator feedback gains. This model is typically used only in determining 
the regulator gains - 5 K- It is not generally used for the design of observer gains (although it 
could be), since that design problem can be reduced to one of lower order, as noted before. 
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Observer Synthesis Model 

The observer uses the measurement vector z (e.g.. measured relative positions and 
accelerations) and the control vector u as inputs to produce observations of the state vector. (In 
general an observer is needed to estimate the unmeasurable system states. Only rarely will all 
system states actually be available for measurement, and never will the disturbance- 

accommodation pseudostates actually be capable of measurement) Kalman-Bucy filter design 

uses knowledge of process noise and sensor noise covariance matrices (and, if necessary, cross- 
correlation matrices) to produce optimal observer gains ( ) L, in the sense of optimality previously 
discussed. If the complete augmented state vector 5 i is to be observed, then an appropriate 
observer synthesis model would be as depicted in Fig. 5. 



Figure 5.— Observer synthesis model, observing all states and pseudostates. 
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(19a) 


S L would be found from the equation 

S L = ( i e ) C r + 5 £,V 2 5 £j X’£.V£j )"' 

where S Q is the unique PD solution to the following ARE, 

5 A 5 Q + s Q 5 A t - = 0 

and where 5 A= 5 A-[ S E } 5 V 2 5 eJ ){ 5 E„ 5 V } 5 eI) C 

! v, = V, - V; 5 £. r ( 5 £, v, s E r ,)-‘ ! e, v/ 

and 5 V It S V 2 , S V 3 , and 5 E, are as defined previously. 

Equations (19b, c, d) reduce to the form 


(19b) 

(19c) 

(19d) 


! /t 'Q+'Q 5 A r -( 5 C ! C r + ’£, S v : 5 £j)(’£. ’V, VM’a } C T + S E, ’v[ s El) 

+ ! £, V, 5 fiJ = O < 19e > 


However, as noted previously, the observer does not need to reconstruct the frequency- 
weighting pseudostates z, and z 2 . This fact will permit an observer of smaller dimension. In this 
case Fig. 6 will be an appropriate observer synthesis model, with pertinent matrices defined as 

follows. 6 L would be found from the equation 

‘L = (‘Q ‘C t +‘E, % ‘El X‘£. % ‘El )"' < 20a > 


where 6 Q is the unique PD solution to the following ARE, 


6 A 6 Q+ 6 Q ° A 1 - °Q 


6ZT 


d c T(. 


6 El 


)-; 6 C 6q + 


i c & J 
C's r 1 


eJ = o. 


(20b) 


and where 

‘A =‘A-(‘E, ‘ V J ‘El){‘E. ‘V,‘El)-‘ ‘C (20C) 

S V, =‘ V, - ‘V 2 ‘ El(‘E , ‘V, ‘ El)-‘ ‘E, *V/ (20d) 

6y (20c) 

6 V 2 = 5 V 2 

^V 3 — *Vj ( 20 ?) 

B E,D, O' 

6 E, = O B, O (20h) 

o o 
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= 


(20i) 


A E,C, O ' 

O A, O 

O O A n 

6 C = [C O E h C h ] GW 


Equations (20b,c,d) reduce to the form 



Figure 6. — Reduced observer synthesis model. 
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Figure 7. — Nominal analysis mode!. 


Nominal Analysis Model 

Once the active vibration-isolation controller has been designed, it can usually be reduced 

in siie by modal truncation and/or balance- and- truncate ("Moore's method, [22]). Then the 

closed loop system can be evaluated, with the controller applied to the actual plant 

Letting the state space system {a fib ,B f , b ,C fib ,D p/b } represent the feedback controller, 

a nominal analysis model can be portrayed as in Fig. 7. Checks on nominal stability can be made 
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by simple eigenvalue checks of this closed-loop system. Nominal performance can also be 
conducted readily using this modeL With the appropriate placement of complex A- blocks to 
represent system uncertainties, one can also evaluate the system for robust stability and robust 
performance guarantees, using the powerful methods of mu-analysis [23]. 


Design Philosophy for the Microgravity Isolation Problem 

With the synthesis framework now in place, as presented above, the designer must choose 
a reasonable strategy in order to use his H* synthesis tools with skill. He must determine what 
states to use, what frequency-weighting and disturbance-accommodation filters (if any) to 
employ, the relative weightings of the resultant frequency-weighted states and control, and the 
relative weightings of the various noise vectors in his system model. The designer must also 
decide which measurements to make, whether to use the full- or the reduced observer synthesis 
model, and whether to conduct the regulator- and observer-gain sub-problems in sequence or in 
combination. 

One primary goal in controller design should be simplicity. The construction of the 
controller will be easiest, and its speed of operation fastest, if its complexity (i.e., the number of 
controller states) is kept to a minimum. To accomplish this aim, the authors recommend that the 
engineer seek to design the controller by starting with basic (unextended) H 2 synthesis and adding 
complexity "one layer at a time." For example, he might first determine if basic LQG is adequate, 
and then add appropriate frequency weighting and disturbance accommodation step by step, 
evaluating after each addition whether or not the design is acceptable. If not, the next layer of 
complexity could be added based on the present design inad e quacies. Once an adequate design 
has been found, it is recommended that the controller order (i.e., number of states) be reduced by 
using modal reduction and /or balance-and-truncate. This "step-up, step-down" philosophy 

should keep controller complexity to a minimum. 

A second fundamental goal should be intuitiveness. Unless the problem is posed in such a 
way as to employ the designer's intuition, he will find it very difficult, especially with a three- 
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dimensional problem, to proceed with any degree of speed. The single most important step 
toward an intuitive problem is the proper choice of plant states. For the microgravity vibration 
isolation problem, the authors believe that a reasonably physical choice is payload relative 
position, payload relative velocity, and payload acceleration. A heavier weighting on payload 
relative position (in the cost functional matrix W ; ), for example, signals the H 2 machinery to 
attempt to increase system stiffness. Similar analogues exist few the other two suggested states, as 
noted before. And at least two of these states are readily measurable for microgravity systems. 
Such state choices, then, allow the designer to assign his weightings with a degree of "physical 
feel," so that extended H 2 synthesis becomes more of a craftsman's design tool, rather than a black 

box for use in a time-consuming trial-and-error approach. 

The designer must also decide whether or not to conduct the regulator- and observer-gain 
design problems independently. The well-known "separation principle" guarantees that for a 
perfectly known system the regulator gains K and the observer gains L can be designed 
independently. One approach, then, would be first to design the regulator to meet the design 
goals, and then to design the observer to produce a state-vector estimate that is "accurate 
enough" over a "sufficient" bandwidth. The frequency-weighting and disturbance- 
accommodation extensions, however, affect the state observations in such a manner that "accurate 
enough" and "sufficient" are quite difficult terms to define. The closed-loop system must be 
analyzed as a whole for this purpose. The existence of an observer bandwidth can also be used to 
enhance overall system performance, so that a full-state-feedback system with inadequate 
performance can actually perform quite well when the observer is added. Since the system 
stability- and performance robustness must ultimately be evaluated for the total closed-loop 
system, it is recommended that the entire controller (Le., the observer-plus-regulator) be designed 
as a unit, rather than in parts. 

For the microgravity vibration isolation problem, studies to date indicate that H 2 synthesis 
extensions arc necessary, if one is to produce a practical control [16]. This being so, there are 
certain frequency weightings that are very reasonable choices to use. At very low frequencies, 
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indirect disturbances (i.e., orbiter positional deviations from a perfectly elliptical orbit), will be 
much larger than rattlespace constraints will allow. In the low-frequency region, then, the 
payload relative displacement should be weighted heavily, and the payload acceleration, lightly. 
These weighting choices could reasonably be expected to call for a controller producing unit 
transmissibility between orbiter and payload, at low frequencies. In the intermediate frequency 
range, where payload acceleration is of most concern, that state should be weighted heavily. At 
higher frequencies, where the plant model is not well known, high control weightings and low 

state weightings should be used to call for reduced control. 

Certain disturbance-accommodation filters, as well, will be appropriate for the problem, 
while others will be inadvisable. From a physical perspective, a more massive experiment would 
be less susceptible either to direct or to indirect disturbances. One could expect, then, that an 
input disturbance filter which models a large direct disturbance would call for a controller tending 
to make the system seem more massive (electronically). On the other hand, an indirect 
disturbance alone (i.e., acting through the umbilical) could be anenuated effectively either by a 
greater system effective mass or by a reduced system effective stiffness. The latter means of 
disturbance attenuation is ineffective for direct disturbances. It also tends to reduce the stability 
robustness of the system. Hence, the designer should be wary of having too large an indirect 
disturbance model. 

Output disturbance accommodation and control noise should be included in the system 
model only if necessary. Research to date does not indicate that either is needed for microgravity 
vibration isolator design. Again, the goal is to achieve a satisfactory controller that is as simple as 

possible. 

Observer design involves the numerical solution of an ARE. An ARE involving matrices 
of smaller dimension will be less susceptible to the numerical difficulties which sometimes attend 
such solution procedures. It is preferrable, then, to use die reduced rather than the full observer 
synthesis model. 
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Concluding Remarks 

Active vibration isolation of microgravity science experiments is a three-dimensional, 

MIMO design problem requiring sophisticated design- and analysis tools. Modem control 
methods provide the most natural setting for handling this problem; and with a suitable choice of 
states, modem-control design can be conducted in a relatively intuitive fashion. The H 2 synthesis 
approach can be extended, using frequency weighting and disturbance accommodation 
techniques, to give the designer great flexibility in building a suitable controller. Implementation 
of these extensions involves a straightforward augmentation of various system matrices, so that 
the ARE-based solution methods of LQG synthesis can be readily applied. Extended H 2 synthesis 
provides the necessary tools for the design of a robust isolation system. This paper has provided a 
general framework for using extended H 2 synthesis to design the controller for such a system. 

In addition to the basic plant model, there are three complementary system models that are 
of use in conceptualizing the synthesis problem. Observer synthesis requires fewer pseudostates 
than regulator synthesis, so two respective system models arc needed to reflect this difference. 

The controller model is developed by combining the observer and regulator models, followed by 
reduction of the controller dimensionality. Attachment of this controller to the basic plant model 
produces an analysis model that can be used, with mu-analysis methods, to evaluate the closed- 

loop system in terms of its stability- and performance robustness. 

This paper has also suggested a general design philosophy, for applying the extended Hj 
synthesis machinery to the particular design problem at hand. In addition to an overall design 
strategy, reasonable state choices were suggested, and basic practical guidelines were given for 
the effective use of frequency-weighting and disturbance-accommodation techniques. 
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